Bifurcation of gap solitons through catastrophe theory 
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In the theory of optical gap solitons, slowly-moving finite-amplitude Lorentzian solutions are 
found to mediate the transition from bright to coexistent dark-antidark solitary wave pairs when 
the laser frequency is detuned out of the proper edge of a dynamical photonic bandgap. Catastrophe 
theory is applied to give a geometrical description of this strongly asymmetrical 'morphing' process. 
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The confinement of optical radiation in periodic media (gratings) with nonlinear response occurs in the form of 
gap solitons (GS), or more properly solitary waves, as first predicted by Chen and Mills Q, and studied extensively 
afterwards J2|-|£|. In Kerr media, th e pro totype model for GS is the following system of hyperbolic PDEs with 
Hamiltonian (conservative) structure P,|10[, which couples the forward u+(z, f) and backward U-(z, t) propagating 
envelopes at Bragg carrier frequency lub 

t(d t + d z )u + + u_ + (X \u_\ 2 + S \u + \ 2 )u + = Q, 

g: a) 

i{d t -d z )u- + u + + {x \u+\ 2 + S |u_| 2 )m_ = 0. 

' Equations (1) have been conveniently written in usual dimensionless units z = TZ and t = TVbT, where Z and T 
are the real-world propagation distance and time, V is the Bragg coupling coefficient, and Vb is the group-velocity 
at Bragg frequency. Moreover S and X are coefficients which specify the relative weight of self- and cross-phase 
modulation, and u± are proportional to real-world envelope amplitudes. Equations (1) are usually analyzed with 
S = 1 and X — 2 (one of the two coefficients can be always set to have unitary modulus by a suitable rescaling of 
the field amplitudes) which describes scalar mode coupling, e.g., in optical fiber gratings [0-p|. Conversely we find 
convenient to leave them as generic coefficients in order to describe both the cases of focusing (5, X > ) and 
defocusing (S, X < |l|) nonlinearity, as well as the two limit cases X = and S — 0, which arise, e.g., when the 
cubic nonlinearity originates from cascading in quadratic media | pd] , ^2[ . Cascading adds improved flexibility since it 
permits to control the sign of the effective Kerr nonlinearity by tuning the wavevector mismatch. We also emphasize 
that, in the case S = 0, Eqs. (1) reduce to the integrable (by means of the inverse scattering method) massive Thirring 
model, and hence the localized waves are strictly speaking solitons. 

Despite the fact that GS arise in a variety of physical settings and models, the importance of Eqs. (1) is threefold: 
(i) they describe with reasonable accuracy optical GS experimentally investigated in fiber Bragg gratings and 
£>Y in corrugated GaAs waveguides S; (ii) though not integrable by inverse scattering method (except for the massive 
,,-C! Thirring case), the model allows to construct the whole family of solitary waves; (iii) Eqs. (1) have allowed to assess 
Oh- the occurrence of peculiar effects such as the onset of oscillatory instabilities ||,[l0) , ultimately related with the absence 
of material dispersion (i.e., second-order derivatives) which distinguishes Eqs. (1) from other soliton-bearing dispersive 
models (e.g., those of the nonlinear Schrodinger type). GS of Eqs. (1) have been studied for more than a decade, and 
both bright and dark Q solutions were reported. Yet, the existence of such GS solutions and their bifurcations 
(how they change qualitatively against changes of parameters) were never investigated to full extent. Here we unveil 
the bifurcation structure of GS, restricting ourselves to subluminal solutions for physical reasons. We show that 
moving Lorenztian GS mark the transition between in-gap bright GS and dark-antidark GS pairs which coexist either 
below or above (depending on the focusing or defocusing nature of the nonlinearity, respectively) the edge of a suitably 
defined dynamical gap. 

Following the notation of Ref. jnj, the entire family of solitary waves of Eqs. (1) can be characterized by seeking 
solutions in the following form 



i+{z,t) = U + ^i{C) exp{^A< + i \f3C + <MC)]} 

(2) 

*_ (z, t) = U. exp{-*Ai + i [/?C + 0- (C)]} 
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where (3 = -/vA plays the role of GS propagation constant, and the intensity 77 and the chirp (nonlinear phase) <fi± 
profiles depend on £ = 7(2 — vt), with 7 = (1 — v 2 )^ 1 ^ 2 being the (subluminal, |i>| < 1) Lorentz factor. Furthermore 

U+ = A^Jj^j and U- = —sA^f^^ account for the velocity induced asymmetry between the foward and backward 

modes, whereas s = sign [X(l — v 2 ) + .5(1 + v 2 )\ is the sign of the nonlinearity which appears only through the 
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common overall coefficient A = 7 _1 \2X (l — v 2 ) +2S (l + v 2 ) . Importantly the entire family of GS is charac- 
terized by two independent parameters, namely the normalized detuning A and soliton velocity v (corresponding, in 
real- world units, to 5lo = TVbA and V = Vbv), which are related to rotational and translational group symmetries 
of Eqs. (1), respectively JlOj . Note that A = corresponds to the Bragg frequency, and v — yields still GS, i.e. 
confinement of light with zero velocity in the lab frame. 

By direct substitution of Eqs. (||) into Eqs. (1), it is readily verified that the intensity 77 and overall phase 9 = <$>+ — <$>- 
obey the following one-dimensional (thus integrable) Hamiltonian system of equations 



• „ • „ dH 
r) = 2r)sm8 = ~-gg, 



= 2S + 2 cos 6-t] = — , 



(3) 

drj 



where the dot stands for d/d(, H = H(r], 9) — 2i]cos6 + 2br\ — rj 2 /2 is the reduced conserved (H = 0) Hamiltonian, 
which now depends on the single parameter 5 = 7A. In Eq. (1) and hereafter we implicitly assume to deal with 
the self-defocusing nonlinearity s = — 1 originally considered in Ref. Q, including the case of vanishing S or X. 
The results can be readily extended to the self-focusing case s = 1 with the substitution 9 — > 2ir — 9 and 6 — ► —S. 
Importantly, the latter condition means that the role of frequencies below and above Bragg frequency must be simply 
interchanged. 

The reduced system (^) permits to find the solitary waves of Eqs. (1) as the separatrix trajectories which are 
homoclinic to (i.e., emanate from and return to) the unstable fixed points r) s ,9 s of Eqs. (^). Given the constraint 
T] > in Eq. (g), these are easily found to be of two kinds: (i) (^cosf^) = (0, —5) for 6 2 < 1, which is associated 
with solutions of the bright type since r)(£ = ±00) = ri s = 0; (ii) (r/ s ,9 s ) = (2 (6 — 1) , 7r) for 5 > 1, which describe GS 
with non- vanishing background or pedestal 7y(C = ±00) = rj s = 2 (5 — 1). For any fixed value of the parameter S, the 
solitary waves correspond, in both cases, to level curves of the Hamiltonian H s = H(r) s , 9 S ). 

Let us first clarify the relation between the existence domain of these two families of GS and the stopband or 
forbidden gap of frequencies exhibited by the grating in the linear operation regime. Bright solutions exist for \8\ < 1, 
which is mapped into the inner domain A 2 + v 2 < 1 of the parameter plane A, v (such representation was already 
adopted in Refs. J10| , ^3| ). This unitary circle can be regarded as a dynamical gap, since in this domain the linear 
problem (i.e., Eqs. (ljwith S,X = 0) yields exponentially damped traveling-wave solutions (u^ oc exp(— iAt + iQ(), 
with Q 2 < 0). On the other hand, solutions with nonzero pedestal exist only outside this dynamical gap or unitary 
circle. 

To clarify further the role of the soliton velocity v, it is important to link the dynamical gap A 2 + v 2 < 1 with the 
bandgap (henceforth termed rest gap to rule out any possible source of misunderstanding) of the stationary linear 
coupled-mode problem. Such rest gap is well known to be related to the reflectivity bandwidth of the grating, in turn 
measurable by means of a tunable cw laser in the laboratory (or rest) frame (z,t). Quantitatively, the rest gap is 
given by the frequency range |<5i| < 1 (\Swi\ < TVb in real-world units) where solutions u±(z,t) = u±(z) exp(— i5\t) 
with frequency detuning 8\ (real-world frequency + 5lj\) are exponentially damped in space, at variance with 
the out-gap case \8\\ > 1 where they become oscillatory jl4). In the velocity- frequency plane (v,S\) the rest gap is 
the (^-independent) lighter shaded domain shown in Fig. 1(a), which has a clear one-to-one correspondence with the 
reflectivity bandwidth reported for comparison in Fig. 1(b). To map the dynamical bandgap in the same parameter 
plane, we need to know what is the actual normalized GS frequency detuning in the rest frame, which is readily 
found to be Si = 7 2 A by grouping in Eqs. (2) phase terms proportional to t. In other words the excitation of a 
GS characterized by the parameters v and A requires a source with detuning = 7 2 A from Bragg frequency. As a 
consequence the dynamical gap where bright GS exist, can be mapped onto the whole shaded domain S 2 < 7 4 (1 — v 2 ) 
or equivalently \S±\ < 7 of Fig. 1. It is clear that, for any given velocity v, this entails a frequency range which is wider 
than the rest gap \S±\ < 1 and reduces to it only in the v = limit. As a consequence bright GS exist for frequencies 
ranging from the rest gap (reflectivity bandwidth) in the v — case, to the whole frequency axis as \v\ — > 1 (i.e., as 
the soliton velocity V approaches the linear group- velocity Vb in the forward or backward direction). 
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FIG. 1. (a) The dynamical photonic bandgap (whole shaded domain) \5i\ < 7 in the parameter plane of velocity and rest 
detuning (v,5i). Such domain is mapped back onto the inner domain bounded by the circle A + v = 1 in the plane (A, v), 
see Fig. 2. The rest bandgap \5\\ < 1 is the smaller region between the two dashed lines 5i = ±1 (light shaded area), and 
corresponds to the bandwidth of the linear reflectivity curve shown with the same vertical scale in (b) for a grating of normalized 
length z L = FL = 4. 



As far as the terminology is concerned, a last important comment is in order. Bright solitary solutions of Eqs. (1) 
are usually classified as gap solitons or Bragg solitons, depending on their detuning 5i being inside < 1, light 
shaded domain in Fig. 1) or outside (1 < |<5i| < 7, dark shaded domain in Fig. 1) the rest gap, respectively. Though 
such a distinction can be useful to locate the operating frequency with respect to the reflectivity bandwidth of the 
grating, it appears otherwise rather arbitrary. First, for any given velocity v there are no qualitative changes of the 
solutions by crossing the boundary \5\\ = 1 between these two regions. Second and more important, the picture of 
Fig. 1 suggests that all the existing bright localized waves are in fact gap solitons if referred to the dynamical gap. 
In other words the effective frequency gap seen by a soliton which moves at velocity v is wider (the faster the soliton 
the wider the gap) than the rest gap measured through the reflectivity. This is also supported by the fact that 8 
can be interpreted as the frequency detuning of the soliton in Lorentz transformed variables £, t with r = j(t — vz) 
[p~5tl - Therefore the rest gap condition \Si\ < 1 must be replaced in moving soliton coordinates by the condition 
\S\ = \S± I/7 < 1 which coincides with the dynamical gap. 

Equations ([}]) are equivalent to the motion of an ideal particle of unitary mass and total energy E in a quartic 
potential well U(rj), ruled by the equation i) = —dU/drj. The kinetic energy is easily obtained in the standard form 
from the first of Eqs. ([}]) by eliminating sin (9 through H, which yields 



77 = V2[E-U(r,)}, (4) 

where U = -25Hr] - \ (4 - AS 2 - H) rj 2 - 8rf + ±?7 4 , and E = -H 2 /2. The GS solutions can be worked out explicitly 
by inverting the quadrature integral obtained from Eq. (Q) with the energy E = —H 2 /2 pertaining to the unstable 
fixed point. We found the following expressions for the case (i) \5\ < 1, entailing bright GS 

4(l-<5 2 ) 

VB ~ cosh(2vT^C) _ ~6' (5) 



tan 



y r T~S 2 sinh(2VT~~^ Q 
-1 + 5 cosh(2VT~S 2 C) 



(6) 



Close to the low- frequency edge of the dynamical gap S ~ — 1, the following approximation of the intensity profile in 
Eq. (§) holds valid [exploit cosh(2a;) = 2 cosh 2 (x) - 1 and 1 - S 2 - 2(1 + 6)] 



r ]B =A{l + S) sech 2 y/2(l + S)C , (7) 



which is characteristic of the one-soliton solution of the focusing nonlinear Schrodinger equation, which provides a 
reasonable description of GS in this region of the gap ||[l6| . 

In the case (ii), i.e. for S > 1 two solutions coexist, being associated to two branches of a double-loop separatrix. 
The first one represents a dark soliton 

x cosh (2VS~1() -VS 
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Odk — tan 1 



4 V / ^T(5 + 1) sinh(2 v / ^ r TC) 



3S 2 + 2Si ~ 3 + 4(5 - 1)<5 cosh(2vtf^T0 + (<5 2 - 25i 
whereas the second one is a bright on pedestal or so-called antidark solution |H 

%/5cosh(2\/^ _ TC) + 1 



1) cosh(4 v / ^TC) 



VAK = 2 (5 - 1 



fl^x = tan" 



\/5cosh (2v / ^ r TC) - l' 

4 v / 5 Tr lV5 sinh(2V^TC) 
2-3S + 6 cosh(4V^TC)" 



(9) 

(10) 
(11) 



These dark and antidark solutions specialized to the zero-velocity case (v = 0) have important implications in terms 
of the stationary (d t = 0) response of the grating, inducing limiting or frustrated bistability, as discussed in Ref. [ pj| . 

Right on the high-frequency edge of the dynamical gap, i.e. for 5 — ► 1, the dark GS vanishes (i~)dk ™> 0), while the 
bright and the antidark solutions have the following common limit 



1 + 



hz = tan 1 



4( 



(2C)' - 1 



(12) 
(13) 



which represents a finite-amplitude moving Lorentzian soliton, i.e., a GS with Lorentzian intensity profile. The 
existence of an exact solution of Eqs. (1) with non-exponcntially decaying tails can be understood from the dynamical 
system (||) as being associated with a degenerate fixed point at the origin with zero eigenvalues. The fact that the 
Lorentzian shape approximates well bright GS close to the upper bound of the gap was noticed earlier under strictly 
stationary conditions (v = 0) ||. 

Viceversa, for S — * —1, i.e., close to the low- frequency edge of the dynamical gap, the intensity of the bright GS in 
Eq. (H) [or, consistently, of its nonlinear Schrodinger approximation (M)] reduces to the following Lorentzian profile 



4(1 + 6) 



4(1 + 5) 



l + 4i±|C 2 1 + 2(1 + <5)C 2 



(14) 



In this case, however, as the stopband edge is approached the Lorentzian GS ( |i4| ) becomes broader and smaller, and 
eventually vanishes identically (t]b — * 0) for S = — 1. The difference between the finite- amplitude [Eq.([T^)1 and the 
vanishing-amplitude [Eq.(p))] Lorentzian GS accounts for an intrinsic asymmetry of the nonlinear grating operation 
with respect to interchange of frequencies below and above the Bragg frequency, respectively. 
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FIG. 2. Existence diagram (defocusing case) for subluminal GS in the region of detuning-velocity (A, v) plane bounded by 
the thin lines at v — ±1. The insets show GS intensity profiles |u±| 2 sampled at the nearest marker (filled circle). Bright 
GS of low-amplitude (LA), high-amplitude (HA), and asymmetric (AS) types exist inside the dynamical gap (unitary circle 
A 2 + v 2 = 1). HA solitons become finite-amplitude Lorentzian (LZ) over the right semi-circle, and then bifurcate into dark 
(DK) and antidark (AK) pairs which coexist for frequencies above the upper edge of the dynamical bandgap. In the focusing 
case the same picture holds with A — » —A. 
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Figure H summarizes the nature of the different GS solutions by reporting the intensity profiles |u±| 2 of qualitatively 
different GS solutions sampled in the parameter plane (A, v). First, it must be noticed that still (v — 0) GS have 
equal intensities |w+| 2 = |w_| 2 as a consequence of the fact that the net photon flux is zero, for both in-gap (bright) 
and out-gap (dark-antidark) solutions. Moreover, let us recall that the amplitude of the in-gap (bright) GS increases 
and their width decreases by spanning the gap from left to right. The symmetry |it+| 2 = |u_| 2 is broken for moving 
GS, which have a stronger component in the direction of motion (i.e., \u+\ > for v > and \u—\ > for v < 0). 
Importantly, bright GS which have high-amplitude (see HA inset in Fig. ||) close to the upper edge of the dynamical 
gap, become Lorentzian GS (LZ inset in Fig. ^]) over the edge A + v 2 — 1, and then bifurcate into dark-antidark pairs 
(DK and AK insets in Fig. |2|) outside the dynamical bandgap. Viceversa, as explained above, the low-amplitude (LA 
inset) GS which exist close to the lower edge of the gap vanishes in the limit 6 — > — 1, and no solutions exist below the 
bottom of the stopband (i.e., outside the circle for A < 0). For a focusing nonlinearity (s = 1) an identical picture 
with A — * —A holds true, meaning that the dark-antidark pairs originate always from the high-amplitude bright GS, 
though, in this case, they now exist below the low-frequency edge of the dynamical gap. 

The bifurcation of GS can be effectively explained in term of the catastrophe theory ]l9|] , and the underlying 
classification of the singularities of smooth functions. The quartic potential U (rj) belongs to the so called cusp 
family (l9) . It can be recast in the following standard form by means of the change of variable r\ — \/2x 4- 25, thus 
obtaining 

U(x) = ix 4 + ^x 2 + bx, (15) 

where a = \/2 (-4 - 2S 2 + H) and b = -8\/2S. A necessary condition for the solitary waves to exist is that the 

potential U(x) has three critical points (i.e., such that dU/dx = 0), thus being of a double- well type. In the 
parameter plane (a, b) this occurs in a domain bounded by the following curve, so-called "bifurcation set" |lS| ], 




which is reported as a thin solid line in Fig. |3|. As shown this curve has the characteristic shape of a cusp. In the 
spirit of the catastrophe theory, we also report in Fig. [|the so-called control line, i.e., how the parameters (a, 6), and 
as a consequence the potential U , vary by changing the single control parameter <5 from large negative values to large 
positive values (indicated in Fig. || by the limit 5 — — oo and d = oo, respectively). In doing so we calculate a with 
H = H s , i.e. the value of the Hamiltonian pertaining to the solitons. 
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FIG. 3. Plot of the bifurcation set (thin line) from Eq. (jl6|) and the control line (thick line) which gives the parametric 
dependence of the coefficients a and b of the potential U(x) in Eq. dli| ) on the control parameter S, when this is varied from 
large negative values (S — —oo) to large positive values (S — oo). The dashed portion of the control line correspond to 
unphysical solutions (77 < 0). Two catastrophes occur at the points 8 = ±1. 

According to our analysis, by varying 5 in the range (—00, 00) the control line remains always inside the cusp, 
indicating the possibility to have solitary waves. Dramatic qualitative changes of the solutions must be expected at 
the catastrophe points 5 = ±1 where the control line becomes tangent to the bifurcation set. In spite of the apparent 
symmetry of Fig. |[ the two catastrophes are the signature of the strongly asymmetrical behavior of GS against 
interchange of frequencies below (5 < 0) and above (S > 0) the Bragg frequency. Indeed, the <5 = 1 catastrophe 
marks the point where the maximum of the potential ( |l5| ) moves from a finite positive value for S > 1 (which makes 
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accessible two distinct asymptotic evolutions inside the two wells, in turn corresponding to the dark and antidark GS), 
to the origin for — 1 < S < 1 (where only one well is accessible for the asymptotic motion towards the origin which 
describes the bright GS). Conversely, the other catastrophe at 6 = — 1 marks the point where these bright GS simply 
cease to exist because the maximum of the potential moves towards negative values. In fact, though the potential still 
has a double- well shape even for — oo < S < — 1, the possibility to have solitary waves is ruled out by the fact that the 
two wells become accessible only with r\ < 0, and hence the solutions are unphysical (recall that r\ is an intensity). 
The dashed line in Fig. || displays that portion of the control line where the solutions are unphysical. 

In summary we have shown that GS solutions of a well-known standard coupled-mode model with Kerr or Kerr- 
equivalent nonlinearity undergo a bifurcation which is strongly asymmetrical with respect to the Bragg frequency. The 
qualitative change of the solutions is explained with a geometrical picture based on the application of the catastrophe 
theory. The bifurcation is marked by the existence of finite-amplitude Lorentzian GS. In this sense it is reminiscent 
of the recently investigated case of localized waves sustained by a gap of full nonlinear origin (2^]. In spite of the 
diversity between the bifurcation discussed here and that of Ref . |2(J] , this suggests that Lorentzian solitons can play 
a universal role in the localization of light in periodic media. 
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